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Coherent population trapping is shown to occur in a driven 
symmetric double-well potential in the strong-field regime. 
The system parameters have been chosen to reproduce the 
CP <-» 3 + transition of the inversion mode of the ammonia 
molecule. For a molecule initially prepared in its lower dou- 
blet we find that, under certain circumstances, the 3 + level 
remains unpopulated, and this occurs in spite of the fact that 
the laser field is resonant with the 0~ <-+ 3 + transition and 
intense enough so as to strongly mix the + and 0~ ground 
states. This counterintuitive result constitutes a coherent 
population trapping phenomenon of nonperturbative origin 
which cannot be accounted for with the usual models. We 
propose an analytic nonperturbative model which accounts 
correctly for the observed phenomenon. 
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Quantum dynamics in symmetric double-well poten- 
tials is important to understand numerous physical and 
chemical processes. A typical example is the tunneling 
dynamics of the hydrogen atoms in the inversion mode of 
the ammonia molecule, which is responsible for the split- 
ting of the vibrational levels [1]. Other examples include 
electron tunneling in quantum semiconductor structures 
[2] or intermolecular proton transfer processes [3] . 

In recent years there has been increasing interest in 
quantum coherence phenomena displayed by atomic and 
molecular systems irradiated with strong laser fields [4]. 
Coherent external fields induce quantum interference ef- 
fects such as coherent population trapping [5], electro- 
magnetically induced transparency [6] or lasing without 
inversion [7]. In particular, in connection with the time 
evolution of a quantum system in a symmetric double- 
well potential it has been shown that under certain cir- 
cumstances an intense laser field can induce coherent tun- 
neling suppression [8]. In this Letter we show that such 
a system can also exhibit coherent population trapping. 
This population trapping phenomenon is nonperturba- 
tive in nature and cannot be accounted for with the usual 
models. We propose an analytically solvable nonpertur- 
bative model which accounts correctly for the essential 
features of the observed phenomenon. 

Specifically, we consider a symmetric quartic double- 
well potential strongly driven by a linearly polarized laser 
field. After appropriate scaling the corresponding dimen- 
sionless Hamiltonian reads 



p2 



(1) 



where the coupling constant A is proportional to the 
laser field amplitude Eo and r = uijj, with wl being the 
laser frequency. The dimensionless parameter a, which 
gives approximately the number of doublets below the 
barrier top, has been chosen to be 1.735. This value 
reproduces, to a good approximation, the effective po- 
tential involved in the inversion mode of the ammonia 
molecule. The laser frequency has been tuned to the 
CP «-> 3 + vibrational transition and its intensity satisfies 
A(0 + \X\ CP) = E /ii2 = 0.357to;l where \x\i is the dipole 
matrix element between the ground states |0 + ) and |CP), 
and all quantities are assumed to be dimensionless. 

Transitions + <-> CP and CP <-» 3 + are dipole allowed 
whereas the + <-> 3 + transition is forbidden. Therefore, 
in the weak-field regime (Eo/iis/wL,, Ao/^l <c 1 with 
Ao being the energy splitting of the lower doublet) and 
for a laser field tuned to the CP «-> 3 + transition, one 
expects the upper level to be populated or not depending 
on whether the molecule is initially prepared in the |CP) 
or |0 + ) state, respectively. The laser intensity considered 
above, however, corresponds to the strong-field regime 
(E /ii2 ~ wl). Under these circumstances the two lower 
levels become strongly mixed and the |CP) state becomes 
highly populated. One then would expect the upper level 
to be populated irrespective of the fact that the molecule 
be initially prepared in the |CP) or |0 + ) state. 

Fig. la shows the evolution of the populations for an 
ammonia molecule initially prepared in its ground state 
|0 + ). These results have been obtained numerically by 
direct integration of the Schrodinger equation. We have 
included the 20 lowest-lying levels, which guarantees con- 
vergence. As is apparent from the figure, the upper level 
|3 + ) remains always unpopulated [curve (3)], and this oc- 
curs despite the fact that the |CP) state becomes highly 
populated [curve (2)] and the laser field directly connects 
this latter state with the upper level. This figure also 
shows that the total population in the lower doublet re- 
mains always close to unity [curve (1)]. Thus, under the 
action of the driving field the initial population oscillates 
rapidly between |0 + ) and |CP) while it remains trapped 
in the lower doublet. This counterintuitive result rep- 
resents a coherent population trapping phenomenon of 
nonperturbative nature which cannot be accounted for 
with the usual models. 

The case of an ammonia molecule prepared initially 
in |CP) is considered in Fig. lb. This figure shows the 
time evolution of both the population of the upper level 
|3 + ) [curve (3)] and the total population of the lower 
doublet [curve (1)]. As before the initial population os- 
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cillates very rapidly between |0~) and |0 + ) (not shown 
for clarity). Now, however, a periodic population transfer 
between the lower doublet and the upper level takes place 
on a different timescale. In fact, apart from the rapid os- 
cillations of the upper level population (which originates 
from population transfers to levels adjacent to |3 + ), as a 
detailed numerical analysis reveals) the behavior of the 
system in this nonpcrturbative regime resembles that of 
the corresponding weak-field regime. 
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FIG. 1. Populations vs. r = uj-^t for an ammonia molecule 
initially prepared in: (a) the |0 + ) state; and (b) the |0~) state. 
Curves (1) give the total population of the lower doublet; 
curve (2) gives the population of the |0~) state; and curves 
(3) give the population of the upper level |3 + ). 

The above results pose two intriguing questions: i) 
Why does the population become trapped when the 
molecule is initially prepared in its ground state? and 
ii) Why does the system behave essentially in a similar 
way both in the nonperturbative strong-field regime and 
in the weak-field regime? In what follows, we propose 
an analytic nonpcrturbative three-level model which can 
give an answer to these questions. 

The most directly involved states, |0 + ), |0~) and |3 + ), 
will be denoted |1), |2) and |3), respectively. The energy 
splitting of the lower doublet is A , and w 3 denotes the 
energy of the upper level. The system Hamiltonian is 



A 

— ((T 2 2 - CTll ) + - ^12 cos ( T ) (0"12 + C21, 



^23 COS (t) (er 2 3 + O32) : 



(2) 



where h = 1; cry = and fiy = with being 

the dipole matrix elements between |i) and \ 

The most rapidly oscillating terms can be absorbed by 
performing the unitary transformation 



i— — (<Ti2 + cr 2 i) sin (r) + i<r 33 T 



U (t) = cxp 
which leads to the transformed Hamiltonian 



(3) 



H' = (A /2) {cos [20(r)] (a 22 - <th) 

+ i sin [20(r)] (cr 2 i - V12)} + (w 3 - ^l)c 3 3 

- n 23 cos (r) { e~ tT (cos [4>(t)] a 23 

- ism [0(r)] CT13) + h.c.} , 
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with (j)(r) = (O12/WL) sin (r). Next, wc expand the time 
dependent coefficients of H' in Fourier series, which al- 
lows us to separate the Hamiltonian into a dominant 
constant contribution Hq and a time-dependent part 
AH'(t). Then, substitution of H' into the evolution op- 
erator of the system shows that when the driving field 
is quasiresonant with the |2) <-> |3) transition and both 
the energy difference Ao and the Rabi frequency fi 2 3 are 
small in comparison with the laser frequency, AiJ'(r) 
becomes a small, rapidly oscillating perturbation which 
can be safely neglected. More generally, it can be shown 
that in the strong-field regime (fii2/wL > 1) and for 
a quasiresonant laser field, AH'(t) becomes negligible 
whenever Ao/^l, ^23 /^l <C ^/£!i 2 /wl- (In our case, 
A q /lu l = 3.28 x Kr 4 ,ft 23 /w L = 0.23, and Q 12 /ujl = 
1.10.) Under these circumstances, the dynamical evolu- 
tion of the system is governed by the Hamiltonian 

H' = ((T 2 2 — On) + (UJ3 - W L )(T33 1^ ((T23 + (T32) , 

(5) 



where the rcnormalized energy difference A and 
Rabi frequency O^g are field-dependent quantities 
defined as A^ = Ao Jo (2£1i 2 /^l) and = 
2^1,(^23/^12) J\ (^12/wl), with J„ being the n th -order 
Bessel function. The Schrodinger equation associated 
with the above Hamiltonian can be readily solved an- 
alytically, and after transforming back one obtains the 
following nonperturbative general solution 



|*(r)>=C 1 (r)|l) + C 2 (r)|2>+C 3 (r)|3> (6) 



where 



d (r) = C[ (r) cos 0(r) + iC 2 (r) sin <j>{r) (7a) 



C 2 (t) = C' 2 (t) cos 4>{t) + iC[ (t) sin 0(t) (7b) 
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C 3 (T)=C' 3 (r)e-* (7c) 

and the C[{t), which are the probability amplitudes as- 
sociated with the Hamiltonian (5), are given by 

C[(t) ^ C' 1 (0)e l ^ T (8a) 




where <5 R = ui 3 — A R /2 — ui L is the renormalized de- 
tuning and £1 R = \J (fi^j) 2 + (<5 R ) 2 is the renormalized 
generalized-Rabi-frequcncy. The physical content of the 
above solution becomes more transparent by considering 
the extended Hilbcrt space of T-periodic state vectors [9] . 
In fact, the basis {|i'(r))} with \i'(t)) = U+(r)\i) turns 
out to be the natural basis to express |W(t)) 

|*(r)) - C[{t)\1'{t)) +C' 2 {t)\2'{t)) + C$(t)|3'(t)> (9) 

As this expression reflects, the dynamical evolution of 
the probability amplitudes corresponding to the renor- 
malized {|«'(t))} states is governed by the Hamiltonian 
H 1 of Eq. (5). Such Hamiltonian has the same form as 
the original Hamiltonian (2) in the limit Q,\2 — > (in 
the rotating wave approximation and in the frame rotat- 
ing with the laser frequency). Therefore, the theory is 
renormalizable in the sense that when analyzed in terms 
of the {|«'(t)}} states, the nonperturbative effects of the 
radiation field on the dynamical evolution of the system 
can be absorbed into the renormalized splitting A R and 
Rabi frequency f2 R 3 , in such a way that the system evolves 
obeying the same Hamiltonian as that of the weak-field 
regime in the rotating wave approximation. In fact, the 
general solution (9) is valid both in the (perturbative) 
weak-field regime (O^/wl, Ao/wl "C 1) and in the (non- 
perturbative) strong-field regime (f^/^L ^ 1)- 

As Eq. (8a) shows, under the action of the coher- 
ent external field, the |l'(r)) state decouples and all the 
population that is initially in state |1) becomes trapped 
in |1'(t)). For a system prepared initially in its ground 
state this implies, in particular, that the upper level |3) 
remains always unpopulated. This occurs in spite of the 
fact that the initial population oscillates very rapidly be- 
tween the |1) and \2) levels, and the latter is directly 
coupled to 1 3) via a laser field tuned to the 2^3 transi- 
tion. This is a coherent population trapping phenomenon 
of nonperturbative origin. 
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FIG. 2. Theoretical predictions for the same situations con- 
sidered in Fig. 1. For comparison purposes, along with the 
analytical results (solid lines) the corresponding exact numer- 
ical results have been plotted again (dotted lines) . 

On the other hand, when the molecule is prepared in 
t = in the state |2), the population difference between 
the upper level and the lower doublet oscillates in time 
as 



W(t) = - 




Fig. 2 shows the theoretical predictions of our model 
for the same situations considered in Fig. 1. Fig. 2a 
corresponds to an ammonia molecule prepared initially 
in its ground state. In this case, the populations Pu(t) 
of the molecular states {\i)} are predicted to be 

P22(t) = sin 2 ( — sinr ] , P33(t) = (11) 

and pn(r) = 1 — P2i{t). Thus, the population of the 
lower doublet, which under these circumstances coincides 
with that of the renormalized state |l'(r)), remains al- 
ways equal to one [curve (1)] and the upper level remains 
unpopulated [curve (3)] , in good agreement with the nu- 
merical results (dotted lines). Fig. 2a also compares the 
analytical result for the population of level |2) [curve (2)] 
with the corresponding exact numerical result. 
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On the other hand, Fig. 2b shows the evolution of the 
populations for a molecule prepared at r = in state |2). 
In this case the population of the lower doublet coincides 
with that of the renormalized state |2'(r)) and, according 
to our model, oscillates in time as [curve (1)] 



p 23 — i5 R p 23 «S^23 (P22 P33) 



Pii(t) + P22 (t) = cos 2 



S R \ . 2 



2loj 



(12) 



while the population of the upper level behaves as 
P33{t) = 1 - Ph(t) - P22{t) [curve (3)]. These results 
are in good qualitative agreement with the corresponding 
numerical results. The main discrepancy between Figs. 
2b and lb comes from the rapid oscillatory behavior of 
the upper level population. As already mentioned, it can 
be shown that this discrepancy, which decreases as the 
laser intensity does, originates from population transfers 
to levels adjacent to the 3 + level, which now have a more 
significant contribution. In fact, if the numerical problem 
is restricted to the three levels most directly involved then 
numerical and analytical results become indistinguish- 
able on the scale of the figures. Our three-level model 
already captures the essential features of the system and 
enables us to understand the dominant behavior of the 
populations in terms of a nonperturbative coherent pop- 
ulation trapping phenomenon. The upper level remains 
unpopulated when the molecule is initially prepared in 
its ground state because such configuration corresponds 
to an initial preparation in the trapping state |l'(r)). 

Next, we will analyze the influence of dissipation on 
the coherent population trapping phenomenon previously 
found. Spontaneous emission effects can be conveniently 
incorporated by assuming that the upper level decays 
radiatively into state \2) with an effective spontaneous 
emission rate T. The dynamics of the system is now de- 
scribed in terms of the density operator p(t) which obeys 
the usual master equation (in which we have retained 
nonsecular terms). By performing the unitary transfor- 
mation (3) one obtains a transformed master equation 
for the density operator p'(t) = U{t)p(t)U + (t), which, 
within the range of validity of our model, leads to the 
following equations of motion governing the time evolu- 
tion of populations and coherences: 



Pll = 2 (1 ~ Atl) ^3 



(13a) 



P22 = ^23 (P32 - P23) + 77 (1 + A o) P33 (13b) 



Pl2 ~ Pl2 *^23/ 3 13 



(13c) 



P'l3 = * (S R + A?) P ' 13 - inl\p' 12 - I (Vl 3 - 7j A 2/>3l) 



P23 



-, A 2P32 



(13e) 



with A„ = J„ (2ft 12 /w L ) (n = 0,2); pf.. = (i\p'(t)\j) = 
(i'(t)\p(t)\j'(t)); andpj^pg. 

As Eq. (13a) reflects, whenever the Rabi frequency f2i2 
coupling the two lower-lying states is nonzero, the upper 
state 1 3) remains unpopulated in the steady state regard- 
less of the initial preparation. As a consequence, in the 
steady state the molecular population becomes trapped 
in the lower doublet and the fluorescence from level |3) 
vanishes. This behavior, which is in sharp contrast with 
the well-known behavior of the system in the SI12 = 
limit, is typical of systems exhibiting coherent popula- 
tion trapping and has its origin in quantum interferences 
involving the two lower- lying levels [5] . 

It is not hard to see from the above equations that, for 
arbitrary external fields (such that 0), the steady- 

state population of \2'(t)) also vanishes so that all of the 
population becomes trapped in the steady state in |l'(i)) 
irrespective of the initial preparation. 

In conclusion, we have shown that coherent popula- 
tion trapping can occur in the nonperturbative regime 
and have proposed an analytically solvable nonpertur- 
bative three-level model which enables us to understand 
the observed phenomenon. Although we have presented 
results for only one field intensity, essentially the same 
behavior, in good agreement with our analytic model, 
occurs in the parameter range O.lir < f2i 2 /wL < 0.57T. A 
detailed account of the model will be given elsewhere. 
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